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Abstract 

In this article we study the geodesic motion of test particles and light in the Einstein-Maxwell- 
dilaton-axion black hole spacetime. We derive the equations of motion and present their solutions in 
terms of the WeierstraB p-, cr- and ^functions. With the help of parametric diagrams and effective 
potentials we analyze the geodesic motion and give a list of all possible orbit types. 


1 Introduction 

The discoveries in the last decades revealed that approximately 95% of the Universe is made of dark matter 
and dark energy. Until today the true nature of dark matter and dark energy remains an open question. 
The string theory, one of the best candidates to solve the problem of quantum gravity, may provide an 
answer. Scalar fields, which arise in some superstring theories, are promising candidates to explain dark 
energy and also dark matter. The low energy limit of the bosonic sector of the heterotic string theory 
is called Einstein-Maxwell-dilaton-axion (EMDA) gravity. This model is a nontrivial generalization of 
Einstein-Maxwell gravity and contains two scalar fields: the dilaton and the axion, which can be related 
to dark energy and dark matter (see e.g. [1]). The four-dimensional EMDA model is described by the 
action 

S = j(^R - , (1) 

where (j) is the scalar dilaton field and n is the axion field, which is dual to the antisymmetric tensor field 
F[ = * Ak/F Ffii, = —\^/^efiuapF°‘^ is the dual of the electromagnetic tensor. A black hole solution 

of the EMDA field equations resulting from this action was found by Garcia, Galtsov and Kechkin [2] . 

A powerful tool to study black holes is geodesics, since the structure of the orbits of particles and 
light provides information about the black hole. The analytical solutions of the geodesic equations are 
useful to calculate the shadow of a black hole or observables like the periastron shift. Those could later 
be compared to observations to test models and theories. 

The first to find the analytical solutions of the geodesic equations in a black hole spacetime was 
Hagihara [3]. He solved the equations of motion in the Schwarzschild spacetime in terms of the elliptic 
WeierstraB p-function. In the Taub-NUT [3], Reissner-Nordstrdm [S], Myers-Perry (with equal rotation 
parameters) [6] and Kerr-Newman [7] spacetimes, the equations of motion are also of elliptic type and 
were solved in the same way. 

However, if black holes in higher dimensions are considered or if the cosmological constant is included, 
the equations of motion are usually hyperelliptic. In a nutshell the integration of the geodesic equations can 
be traced back to the Jacobi inversion problem and the solutions are given in terms of the tr-function. This 
was done for the four-dimensional Schwarzschild-de Sitter spacetime 019] as well as for higher-dimensional 
Schwarzschild, Schwarzschild-(anti)de Sitter, Reissner-Nordstrom and Reissner-Nordstrom -(anti) de Sitter 
spactime [10]. Likewise the geodesic equations were solved in the Kerr-(anti) de Sitter spacetime [11], in 
the higher-dimensional Myers-Perry spacetime |12] and in the Hofava-Lifshitz spactime m- In five¬ 
dimensional black ring spacetimes, the equations of motion could be solved analytically in special cases 
[ms]- Moreover, the motion of test particles was studied in various black string spacetimes [MEZKIHl 
[Ta|20]|2l] including field theoretical cosmic string spacetimes [nniinniEi] and black holes pierced by a 
black string [25l |26] . 


In this article we study the geodesic motion in the EMDA black hole spacetime. We analyze the 
possible orbit types using effective potential techniques and parametric diagrams. Furthermore we present 
the analytical solutions of the equations of motion for test particles and light. The equations of motion 
are of elliptic type and the solutions are given in terms of the Weierstrafi p-, a- and ^-functions. 


2 The Einstein-Maxwell-Dilaton-Axion black hole 


The metric of the stationary axisymmetric EMDA black hole is given by [2] 


with 
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dp^ , 


E = — 2mr + c? , 

A = — 2/?r + cos^d . 


( 2 ) 


(3) 


In comparison with the original solution of [2], here we chose b = 0. If the parameter b is different from 
zero, the solution carries a NUT charge. 

The parameter /3 is the dilaton charge, m is related to the ADM mass M = to — /3, and a corresponds 
to the angular momentum J = aM. The electric charge of the black hole is Q = \/—2ijjf3M, where 
a; = exp(20o) and (j)o is the asymptotic value of the dilaton field (see [2]). Demanding that Q is real and 
M > 0, we can conclude that /3 < 0. If /3 = 0, the solution reduces to the Kerr spacetime. 

The horizons are given by E = 0 and can be written in terms of the ADM mass 

r±=M + P±^J{M + pf-a^. ( 4 ) 

The condition < (M + (3)'^ has to be fulfilled in order to get real values of r±. In contrast to the Kerr 
spacetime, here it is possible that and r_ are negative. In this article we will only analyze the case 
of positive horizons in detail, but nonetheless the analytical solutions of the geodesic equations given in 
later sections are valid in all cases. 

The ergoregion is determined by E — sin^d = 0, for example in the equatorial plane (d = ^) the 
boundary of the ergoregion is at r = 2 (M + /3). Inside this region the test particles have to rotate in the 
same direction as the black hole. 

The singularity is obtained for A = — 2/3r + a?' cos^d = 0, since the Kretschmann scalar diverges 

for A = 0. Due to the influence of the dilaton parameter /3 the shape of the singularity is different from 
the ring singularity of a Kerr black hole. At r = 0 and d = ^ we still have the ring singularity, but for 

negative r the singularity extends up to r = 2/3 (there is no singularity at positive r). Its shape varies 

from a deformed torus to a complex structure of two closed surfaces depending on the parameters a and 
/3, see figure [1] for a two-dimensional projection. However, only geodesics with r = 0 and d = ^ end in the 
singularity. Geodesics coming from positive r can cross r = 0 if d and continue to negative r. For 
negative r values the geodesics enter a universe where gravity is repulsive (which can be seen in plots of 
the effective potential and the orbits in later sections); therefore orbits with r < 0 cannot hit the singularity. 


The metric is given in Boyer-Lindquist-like coordinates (r, d, p), which can be transformed to Cartesian 
coordinates by 

X = a/( r^ -I- a?) sind cos p 

y = \/ -I- o^) sin d sin p (5) 

z = r cos d . 


dS 1 dS dS 

-h -q^ - 

dr 2 dx’' 


2.1 The geodesic equations 

The Hamilton-Jacobi equation 


( 6 ) 















(a) a — 0.1 and /3 = —0.08 


(b) a — 0.08 and /3 = —0.08 


(c) a = 0.078 and /3 = -0.08 


Figure 1: Plots of the singularity given by A = 0, note that r < 0 in the figures. 


can be solved with an ansatz for the action S 

S = ^St — Et + Lip + Sr(r) + S^{'d ), 


(7) 


where r is an affine parameter along the geodesic, E is the energy and L is the angular momentum of the 
test particle. The parameter 5 represents the mass of the test particle and is equal to 1 for a particle with 
nonzero mass and equal to 0 for light. 

The Hamilton-Jacobi equation (jS]) separates with the help of the Carter m constant K and yields 
four differential equations 
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The polynomial R and the function 0 are given by 
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( 10 ) 

( 11 ) 


( 12 ) 


with 


S = — ^1 + 2 / 3 ^ f + , 

A = — 2pf + d^ . 


(13) 


For convenience we introduced dimensionless quantities by scaling with the ADM mass 
^r~t_T^L_ a~/3^iF 

We also employed the Mino time [28] y as Ady = df. 


(14) 


3 Classification of the geodesics 

The function 0 in equation (|9|) and the polynomial R in equation (|8|) define the characteristics of the 
geodesics. In this section we will therefore study the function 0 and the polynomial R to determine the 
























possible orbits of light and test particles in the EMDA spacetime. In [3D] a similar analysis was done 
for the Kerr spacetime; the Kerr (anti) de Sitter spacetime and the rotating black string spacetime were 
discussed analogously in m and ED, respectively. 


Geodesic motion is possible if 0 > 0 and i? > 0, then real values of the coordinates i) and f are 
obtained. From 0 > 0 it follows immediately that K > 0. It K — 0, the motion is conhned to the 
equatorial plane {d = ^) and L = aE; otherwise an orbit is not possible |5D] . 

As discussed before, a geodesic hits the ring singularity if f = 0 and simultaneously = ^. In this 
case we have 

0 f = K — [aE — L)^ > 0 and 

\ 2 ) ’ - ( 15 ) 

i?(0) = -a^[k - {hE - if] > 0 , 

which yields a condition for orbits ending in the singularity: K = {aE — Lf. This orbit also resides in 
the equatorial plane. Furthermore equation m leads to the following conclusions: 

1. AT > {aE — Lf: The geodesics cross = f but do not cross f = 0. 

2. K < {aE — Lf: The geodesics do not cross ^ but it is possible to cross f = 0. 


3.1 The t?-motion 


The i?-equation is not influenced by the dilaton charge and therefore is the same as in the Kerr black hole 
spacetime. Here we will briefly present the main results of I3Q] and m- 
First we substitute = cost'd (with z/ G [0,1]) into the function 0(d): 


e{iy) = k- Sa^iy - a^E‘^{l - v) + 2dEL 


'k 

(1-z.)- 


(16) 


Since the zeros of 0 are the turning points of the latitudinal motion, we are interested in the number of 
zeros in the interval [0,1]. This number changes if a zero crosses 0 or 1, or if a 0 has a double zero. From 
Q{v = 0) = 0 and Q{v = 1) = 0, which is only possible for L = 0 due to the pole at = 1, we get two 
boundary cases of the d-motion: 


L = aE ± \/k and 


k = 6k A L = o. 


(17) 


To remove the pole of 0 at z^ = 1 we consider 



0'(i/) = k{6 - E‘^)k + [2hE{dLE -L)- 5k -kf + k- {aE - Lf , 

(18) 

where 0(z/) = 

T^0^(z^). Double zeros fulfill the conditions 



0'(z/)=O and 

du 

(19) 

which gives 

- _ E±{k6- k)^/E‘^ -6) 

2dS 

(20) 


Therefore double zeros only occur for > 6. 


Parametric diagrams can be drawn with the help of the equations m and (uni) (see figure [3]). These 
diagrams reveal four regions with different numbers of zeros. In regions (a) and (d) orbits are not pos¬ 
sible since 0(z/) < 0 for all i' G [0,1]. There is one real zero in the interval [0,1] in regions (b), here 
k > {aE — Lf and the geodesics cross d = ^. Two real zeros in the interval [0,1] can be found in region 
(d), where k < {aE — Lf and d = ^ is not crossed. Crossing f = 0 is possible in region (d) but not in 
region (b). 

Consider the special case Z = 0, in which 0 has a single zero at z/q = • For k = 5aL this zero 

is located dX v = 1 and for k = kE^ it is at z^ = 0. For k f Sk, there are two cases 







K > da?: uq G [0,1] if if < and Q{v) > 0 for all v G [0,1] if X > a^E^. 

2. K < 5a?'. vq G [0,1] if if > a?E'^ but Q{v) < 0 for all v G [0,1] if ^ < a?E'^. 

3.2 The f-motion 

The polynomial R determines the possible orbit types, since its zeros are the turning points of the geodesics. 
One can distinguish between bound orbits moving between two turning points and escape orbits, where 
the test particles approach the black hole, turn around at a certain point and then escape to infinity. 
So-called terminating orbits end in the singularity at f = 0 and d = ^, here if = {dE — is required. 
Sometimes the test particles cross both horizons twice or even multiple times. Since these orbits enter 
another universe, they are called two-world or many-world orbits. Similar to the Kerr spacetime, the 
f-coordinate is allowed to take negative values. An orbit that crosses f = 0 is called transit or crossover 
orbit m- Let us first introduce a list of all possible orbits (for r± > 0) before analyzing the f-motion in 
detail. 

1. Transit orbit (TrO) with range f G (— 00 , 00 ). 

2. Escape orbit (EO) with range f G [ri,oo) with ri > f_|_, or with range f G (—oo,ri] with ri < 0. 

3. Two-world escape orbit (TEO) with range [ri,oo) where 0 < ri < r_. 

4. Crossover two-world escape orbit (CTEO) with range [ri,oo) where ri < 0. 

5. Bound orbit (BO) with range f G [ri,r 2 \ with 

(a) ri, r 2 > r+ or 

(b) 0 < ri,r 2 < r_ or 

(c) ri,r 2 < 0 if /3 0. 

6. Many-world bound orbit (MBO) with range f G [ri,r 2 ] where 0 < ri < r_ and ^2 > r_|_. 

7. Terminating orbit (TO) with ranges either f G [0, 00 ) or f G [0, ri] with 

(a) ri > f+ or 

(b) 0 < ri < f_ or 

(c) ri < 0 if ,5 0. 

Note that the bound orbit with ri, r 2 < 0 and the terminating orbit with ri < 0 can only exist in the 
EMDA spacetime with (3^0, but not in the Kerr spacetime. 

The type of an orbit depends on how many real zeros the polynomial R has. The number of zeros 
changes if double zeros appear, that is 


i?(f) = 0 and = 0 . (21) 

dr 

In the parametric L-A^-diagrams based on (EID five regions with different numbers of zeros can be dis¬ 
tinguished. In region (I), the polynomial R has no zeros. Region (II) has one positive and one negative 
zero, whereas region (V) has two positive zeros. In the remaining regions (III) and (IV) four real zeros 
are possible. Depending on the parameters of the metric and the test particle, region (III) has either 
four positive or two positive and two negative zeros. Region (IV) allows the following constellations: one 
negative and three positive zeros, three negative and one positive zero, or even four negative zeros. Eigure 
[5] shows three examples of the parametric L-A^-diagrams of the f-motion combined with the parametric 
diagrams of the d-motion. 

Eurthermore, the possible orbits can be visualized by an effective potential, which consists of two parts 
V+ and V -, determined by the equation 

R{r) = (r^ - 2/3f -f (A - V+) {E - V_) , 


(22) 



0.3 



Figure 2: Combined parametric L-£^^-diagrams of the f-motion (green lines) and i^-motion (blue lines). In region (I) the polynomial 
R has no real zeros. Two zeros are possible in the regions (II) and (V). In the regions (II) and (V) R has two zeros. The 
?9-equation does not allow geodesic motion inside the grey areas (a) and (c). The orbits cross 'd — ^, but not f = 0 in 
regions marked with (b). In contrast, in regions marked with (d) r = 0 can be crossed but -i? = ^ is never crossed. 


and thus 


^± = 




— 2/3f + a? 


(23) 


The area in between the two parts of the potential has R{f) < 0, so there geodesic motion is not possible. 
In the limit f —>■ ±oo, V+ converges to and V- converges to —\/S. If the parameter a or L changes its 
sign, the effective potential is mirrored at the f-axis. and V~ meet at the horizons r± and if if < 6/3'^ 
they also meet at 


r\b = /I± W/32- 


K 


(24) 


leaving a gap between fa and ft. Another peculiarity of the EMDA spacetime is the divergencies in the 
effective potential for at 


fc,d = /3 ± 



(25) 


Taking all the information into consideration, we can now identify all possible orbits in each region of 
the parametric L-A^-diagrams (below it is assumed that < fi+i): 


1. Region (I)(d): There are no real zeros and R{f) > 0 for all f. The only possible orbit is a transit 
orbit which crosses f = 0 , but not f) = ^. 

2. Region (II)(b): R{f) has a negative zero fi and a positive zero f 2 , and R{f) > 0 for f G (—oo,fi] 
and f G [f 2 , 00 ). The geodesics cross the equatorial plane. For f < 0 an EO is possible and for f > 0 
a TEO is possible, whose turning point can coincide with the inner horizon f_. In the special case 
f 2 = 0, the TEO turns into a TO and the motion is confined to the equatorial plane. 

3. Region (II)(d): R{f) has two negative zeros fi, f 2 with R{f) > 0 for f G (—oo,ri] and f G [f 2 ,oo). 
There is an EO and a CTEO; this time the equatorial plane is not crossed. 

4. Region (III)(b): R has four zeros fi, f 2 , rs, r 4 and R{f) > 0 for f G [fi,r 2 ] and f G [fa, 1 = 4 ]. If 
all four zeros are positive, there is a MBO between ri and r 2 (one of the zeros can coincide with a 
horizon), and a BO between fa and f 4 . It is also possible for /3 ^ 0, that fi, f 2 are negative and 
fa, f 4 are positive. Then there is a BO between fi and f 2 , and an MBO between fa and f 4 (one of 
the zeros can coincide with a horizon). In the special case K = {aE — L)^, the BO turns into a TO 
since f 2 = 0. In this region the geodesics cross = §■ 

5. Region (IV)(b): R has four zeros fi, f 2 , fa, f 4 and R{f) > 0 for f G (—oo,fi], f G [f 2 ,fa] and 
f G [r 4 ,oo). Either one zero is negative and three zeros are positive, or three zeros are negative 
and one zero is positive. The latter is only possible if /I ^ 0. EOs exist for f G (—oo,fi] (this one 
resides at negative f) and f G [f 4 , 00 ) if f 4 > f+. If f 4 < f+, then there is a TEO. Depending on the 






















position of the zeros, various orbits are possible between f 2 and r^. There can be MBOs (one of the 
zeros can coincide with a horizon), BOs are hidden behind the horizons with f 2 ,f 3 < f_ and even 
TOs are present ii K = (aE — L)^. For ,5 ^ 0 it is possible that f 2 ,'f 3 < 0 and therefore a BO for 
negative f exists. In this region the geodesics cross = f ■ 

6. Region (IV)(d): This regions only occurs if ,0 ^ 0; here R has four negative zeros fi, f 2 , ta- 
R{r) > 0 for f G (—oo, fi], f G [f 2 , f 3 ] and f G [r^, oo). An EO, a BO and a CTEO exist. Here f = 0 
can be crossed, but not'd = f. 

7. Region (V)(b): R has two positive zeros fi, f 2 and R{r) > 0 for f G The possible orbit is 

a MBO, where one or both (L = 0) zeros can coincide with the horizons. If for K = {aE — L)^ we 
have fi = 0, the orbit turns into a TO. The geodesics cross the equatorial plane. 

Table [T] shows an overview of all possible orbit types in the EMDA spacetime. In figure [3] examples 
for the corresponding energy for each orbit type are plotted together with the effective potentials. Note 
that the repulsive nature of the singularity at negative f can be observed in figure [3Ke) andEDf), where 
potential barriers prevent the test particles from falling into the singularity. The corresponding plot of the 
singularity, which in this case consists of two closed surfaces, is shown in figure [ijc). The energy chosen 
for the orbit type G allows an EO at negative r, a BO at negative f enclosed by the surfaces of singularity 
and a CTEO. 



(d) 5 = 1, a = 0.4, /3 = -0.08, K = 0.006, (e) 5 = 1, a = 0.078, 
i = 0.1 K = 0.008, L = 0.1 


-0.08, (f) 5 = 1, a = 0.078, /3 = -0.08, 

A = 0.006, L = 0.1 


Figure 3: Plots of the effective potentials with examples of energies for the different orbit types of table[2 The blue curves correspond 
to V+('r) and the green curves correspond to V_(f). The grey area between V+ and V_ marks a zone, where geodesic 
motion is not possible since there R{r) < 0. In the dashed area, the motion is forbidden by the i^-equation (©(t 9) < 0). 
The red dashed lines represent the different energies. The red dots are the turning points of the orbits and accordingly 
the zeros of the polynomial R. 
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Table 1: Types of orbits of light and particles in the EMDA spacetime. The range of the orbits is represented by thick lines. The 
turning points are marked by thick dots. The two vertical double lines indicate the position of the horizons and the single 
vertical line corresponds to f = 0. The column S shows whether the orbit is possible for particles (1) and/or light (0). 


4 Solution of the geodesic equation 

In this section the analytical solutions of the geodesic equations ([5 |) - (|TT|l are given in terms of the Weier- 
straB p-, (T- and C-functions. 
















































4.1 The f-equation 

The fourth order polynomial R = Ylt=i right-hand side of equation ([5]) with the coefficients 

ai = E'^ -5 


as = -4/3 (i3^ - (5) -I- (5 

02 = 2d.E {aE-L^-k- - 2(5 (l 2^) -h 4/3^^;^ 

ai=k(l + 2^) -h 25^(5^ - 4h^T; [aE - 


ao = a 



(26) 


can be reduced to order three by substituting f = ±i -|- f/{ (where fa is a zero of R): R' = X)i=o 
Then x = -^ {^y — y) converts R' into the WeierstraB form PJ and ([ 8 ]) acquires the form 


= 4y^ - ff 2 y - ffs = Psiy)^ 
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The elliptic differential equation (l?7l) is solved by the WeierstraB p-function 


7/(7) = P (7-7in; 92.55) , 


(27) 

(28) 
(29) 


dy 


where -C = 7,. + C .»"■ - -w. . 12 

the f-equation ([ 8 ]) 


with yin = ± 4 ^—h Finally, resubstitution yields the solution of 
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(30) 


4.2 The t?-equation 

First we substitute r = cost'd to simplify the i?-equation 
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This transforms the right-hand side of the "d-equation into a third order polynomial X)i=i with the 
coefficients 
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Since equation (1^ is solved by the WeierstraB p-function, the solution of equation ([9]) is: 

i9(j) = arccos (4p(7 - Ti"; 92.93) “ j) 


(32) 

(33) 
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(37) 


where 7 " = jin + /“ 


du 
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4.3 The (y9-equation 

The i^-equation (US can be rewritten using the f- and d-equations (HI and m 


dtp = 


E L 


r — 2f3f + dr\E — dL 


df 


-i— iaE sin^i? — Z) 

in^d V J 


d^ 
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(38) 


Integrating this equation gives an f-dependent integral If and a d-dependent integral which can be 
solved separately 


p-pQ = 
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S 1 


— 2/3f + a?) E — aL 


df 

J y/R 
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sin^i? V 
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(39) 


First, let us consider the integral If and substitute f = ±—+ f/j. A partial fraction decomposition 

4y-ir 

yields 


If = 


' 2/in 




C, 


dy' 




(40) 


where the first order poles of If are pi = and P 2 = ’ which correspond to the 

horizons. The sign of ±63 is determined by the chosen sign in the substitution f = ±—+ rji. Ci are 

constants from the partial fraction decomposition depending on the parameters of the metric and the test 
particle. A further substitution y = p {'j — 92 ^ 93 ) =■ Pf{v) with u = 7 — results in 


If = 


re +v_^ 


) -Pz 


du'. 


(41) 


Analogously the integral can be brought into this specific form with the substitutions if = cos^ d, then 
^ ^ - f) and lastly u = pY - 7 ;" ; g^gt) =: p^(f) with f = 7 - 7 /;^ 


U = 


dE - 


C 3 L 


1 


4 psiv') - q 


df', 


(42) 


where g = ^ + jl is a first order pole. 

If and I^ are elliptic integrals of the third kind and can be solved with the help of the p-, a- and 
C-functions as shown in 1113112. Finally, the solution of the p-equation cni) is 


a 

C 3 L 


pin) = CoY - Win) + ^ YJlYd ( 2 C?(wj)(w - Win) + In 
- dE{v - fin) - 


af{v - Vi) af{v + Vi) 

— In ■ 


‘ip'A^g) 


O-f(win-Wi) Cri=(uin + Wi) 

a^iv-Vq) a^{v + Vq) 


(43) 


2(Yvg)iv - fin) + In ^7: - In 7, ' 

(^■9 (Win -Vq) (Jf) (Uin + Vq) 


Pin 


with p^ = pf (wi), q = p^(fq), w = 7 - 7 i'„, f = 7 - 7 ," and 


pn(w) 

= P(w,d5.d3) . 

Pi9(w) 

= piv,gt,gt), 

Cf(w) 

= C(w,dZd3) > 

Ci9(w) 

= C(w,d7d3) . 

wf(w) 

= cr{v,gY93) > 

a^iv) 

= cr(w,g7d7 • 


(44) 


4.4 The /-equation 

The t-equation m can be rewritten using the f- and d-equations ([ 8 ]) and 


dt = 


Y — 2(31 + d 


2 r 


Y — 2/31 + d^] E — dL 


' df 
- y/R 


— d I dE sin^'d — L 


dd 


(45) 


Integrating this equation gives an f-dependent integral and a d-dependent integral as in the case of the 
p-equation m 


i i'm — 


^ 0 + 
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pf(w') -p* 


du' — 


d'^E - aZ -I- y-^—pi?(f) ) df', 


(46) 




























with the first order poles pi = P'^ ~ ^^i 2 (f ^-Pr)^ ' Here we used the same substitutions 

as in section oi and applied a partial fraction decomposition which gave rise to the constants C'. Again 
we find elliptic integrals of the third kind that can be solved with the help of the p-, a- and ^-functions 
as shown in HHllS] ■Finally, the solution of the t-equation m is 


1, \ n'I \ fo/- f ^ , 1 af{v + Vi) 

t{"f) = Cq{v - Vin) + 2^ — rTTT ( - Uin) + In —;-r - In 


7^1 

f-iTn ~T C2\ 

- ( a^E - aL + —j {v - Wm) - 


C3 


CTf {Vin - Vi) 

iC^iv) - Ci?(^in)) + ^in 7 


with p, = pf{vi), u = 7 - 7 i'„, h = 7 - 7 " and 



Pfi'ij) =p{v,gl,gl), 

P^fiv) 

= Piv,9t,9t), 

(f(v) = Civ, 92, 93 ), 

Ci)iv) 

= Civ, 9 t,gt), 

af{v) =a{v,g^,gg, 

asiv) 

= (^iv,9t,9t) ■ 


a-f (vin + Vi) 


(47) 


(48) 


5 The orbits 

In this section we use the analytical solutions of the equations of motion to plot some examples of test 
particle orbits in the EMDA black hole spacetime, see figure |31-figure [51 A bound orbit and an escape orbit 
are depicted in figure Idja) and figure |4Kb) respectively. A terminating orbit, which lies in the equatorial 
plane, is shown in figure IHJc) . There one can also observe that the orbit (blue line) is discontinuous when 
it crosses a horizon. This is due to divergencies in the p-equation, which occur at the horizons. 

Figure El shows three orbits that cross both horizons, but do not end in the singularity. On the two- 
world escape orbit (figure EKa)) and the many-world bound orbit (E{b)) the test particles emerge into 
another universe each time the horizons are crossed twice. In figure Elb) the effect of the ergosphere is 
visible, since here the angular momentum of the black hole and the angular momentum of the particle 
have opposite signs. If the test particle enters the ergosphere, it is forced to rotate along with the black 
hole and therefore changes its direction. The particle on the bound orbit in figure |4l(a) has also an angular 
momentum with a different sign than the rotation of the black hole, but this motion takes place far away 
from the ergosphere and maintains its direction. 

The test particle on the transit orbit depicted inEJc) crosses f = 0 (but not the equatorial plane) and 
travels into a universe where gravity acts repulsive. Therefore in the orbit plot it looks like the particle is 
reflected at that point. 

In figure E] two bound orbits (figure EKa) and Elb)) and an escape orbit (figure EKc)) with negative f 
around a naked singularity can be seen. The singularity at negative f acquires different shapes depending 
on the parameters of the metric; see also section [2] Bound orbits at f < 0 are always hidden inside the 
“walls” of the singularity. This region enclosed by the surfaces of the singularity can never be reached 
from the outside, nor can a particle leave the inside, since the singularity is repulsive here. This is also 
observed in the plots of the effective potential in figure|3l(e) and|3l(f), which are similar to the corresponding 
potential of the bound orbit shown inEl^b). Although the region inside the singularity can never be reached 
from the outside, we display the orbits for the sake of completeness. 

6 Conclusion 

In this article we studied the EMDA spacetime with the help of geodesics and presented the analytical 
solutions of the equations of motion. The geodesic equations are of elliptic type and can be solved in terms 
of the WeierstraB p-, a- and ^-functions. The dilaton parameter /3 influences the r-equation, p-equation 
and t-equation, but the d-equation is not affected and is exactly the same as in the Kerr spacetime, which 
is recovered by setting ,5 = 0. The dilaton charge changes the shape of the singularity. At f = 0 and 
■d = ^ we still have the ring singularity, but for negative f the singularity forms a toruslike structure or 
two closed surfaces depending on the dilaton charge and the rotation parameter. Similar shapes of the 
singularity of a black hole in Einstein-Maxwell-dilaton-axion theory have been observed in m- Since 
gravity becomes repulsive for f < 0 , geodesics can only hit the singularity at r = 0 and = f (as in the 
Kerr spacetime, these geodesics lie in the equatorial plane). 









(a) 5 = 1, a ^ 0.4, L ^ -0.5, ^ ^ 5, (b) 5 ^ 1, a = 0.4, L ^ 0.8, K = 2, (c) 5 = 1, a ^ 0.4, L = 1, K 0.36, 
0 — —0.08, E — 0.97: Bound orbit 0 — —0.08, E — 1.002: Escape orbit 0 — —0.08, E — 4: Terminating orbit 

Figure 4: Different orbits of test particles in the EMDA black hole spacetime. The blue lines represent the orbits and the ellipsoids 
indicate the positions of the outer horizon r^. and the inner horizon r_. The red circle in figure (c) is the singularity at 
r = 0 and -d — ^. 



(a) 5 = 1, a = 0.4, L — 1, K — 0.5, (b) <5 = 1, a = 
0 — —0.08, E — 5.5: Two-world escape 5, ^ — —0.08, 
orbit bound orbit 


0.4, L = -1, a: = (c) 4 = 1, a = 0.4, L = 0.5, K 
E = 0.8: Many-world 0 — —0.08, E — 5: Transit orbit 


1 , 


Figure 5: Different orbits of test particles in the EMDA black hole spacetime. The blue lines represent the orbits and the ellipsoids 
indicate the positions of the outer horizon r+ and the inner horizon r_. 



(a) 5 = 1, a = 0.2, L = 1, a: = 0.008, (b) 5 = 1, a = 0.078, L = 0.15, AT = 5, (c) 5 = 1, a = 0.35, L = 1, A' = 2, 

0 — —0.08, E — 4.7: Bound orbit with 0 — —0.08, E — 1.1: Bound orbit orbit 0 — —0.08, E — 8.5: Escape orbit with 

r < 0 with r < 0 r < 0 


Figure 6: Different orbits of test particles in the EMDA black hole spacetime. For the orbits depicted here, the coordinate r is 
negative and the singularity is not covered by a horizon. The grey structure is the singularity (compare section[^ figure 
0. In figure (a) and (b), we removed a part of the singularity from the plot, so that the orbit inside becomes visible. 
















To determine the possible orbits in the EMDA spacetime, parametric L-E^-diagrams and effective 
potential techniques were used. We found terminating orbits that end in the singularity and lie in the 
equatorial plane, escape orbits and bound orbits. In this spacetime escape orbits are possible both for 
positive and negative r. Bound orbits occur for f > f+ and for 0 < f < f_ hidden behind the inner 
horizon, but also for f < 0. These bound orbits at negative r do not occur in the Kerr case /3 ^ 0 and 
are hidden inside the “walls” of the singularity. Other possible orbits are two-world escape orbits and 
many-world bound orbits which emerge into other universes every time the horizons are crossed twice. 
Additionally it is possible to cross f = 0 and travel into an antigravity universe, resulting in transit orbits 
with the range —oo < f < oo and crossover two-world escape orbits which cross both the horizons and 
f = 0. 

The exact orbits and their properties can be calculated using the analytical solutions of the geodesic 
equations. Additionally, observables like the periastron shift of a bound orbit or the light deflection of 
an escape orbit can derived with formulas analogous to m- The geodesics of photons can be used to 
calculate the shadow of the black hole. In the EMDA spacetime the shadow was studied numerically in 
[32] . Hopefully one day all this could be compared to observations. 

Eor future work it might be interesting to extend the equations of motion and their solutions to 
electrically and magnetically charged particles. Then the influence of the dilaton on the charged particles 
could be studied and compared to the particle motion in e.g. the Kerr-Newman spacetime. The analytical 
solutions to the geodesic equations of charged particles in the Kerr-Newman spacetime were found in 

[7]. Moreover, the NUT parameter could be included by setting b ^ 0 in the spacetime of [2]. Another 
option is to consider particle motion around a more general black hole spacetime containing more charges 
and a cosmological constant. An interesting candidate for this is the dyonic AdS black hole in maximal 
gauged supergravity found by Chow and Compere [33]. The EMDA black hole studied in the present 
paper represents a special case of the spacetime in |33j . 
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